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Abstract. We briefly introduce hysteresis in spatially extended systems and the 
dynamic phase transition observed as the frequency of the oscillating field increases 
beyond a critical value. Hysteresis and the decay of metastable phases are closely 
related phenomena, and a dynamic phase transition can occur only for field am- 
plitudes, temperatures, and system sizes at which the metastable phase decays 
through nucleation and growth of many droplets. We present preliminary results 
from extensive Monte Carlo simulations of a two-dimensional kinetic Ising model in 
a square-wave oscillating field and estimate critical exponents by finite-size scaling 
techniques adapted from equilibrium critical phenomena. The estimates are consis- 
tent with the universality class of the two-dimensional equilibrium Ising model and 
inconsistent with two-dimensional random percolation. However, we are not aware 
of any theoretical arguments indicating why this should be so. Thus, the question 
of the universality class of this nonequilibrium critical phenomenon remains open. 



1 Introduction 

Hysteresis often occurs when a bistable or multistable system is driven by 
an oscillating force which varies too fast for the system to respond without 
a phase lag. In fact, the word was coined by Ewing from the Greek husterein 
(vaTepeuj) which means "to be behind" ||]. Examples include ferromagnets 
and ferroelectrics in AC fields, electrochemical cyclic-voltammetry experi- 
ments, and nonlinear elastic media under oscillating stress, just to mention a 
few. The most familiar representation of hysteresis is probably the (usually) 
closed curve obtained by plotting the system response versus the applied 
force. Examples of such hysteresis loops are shown in Fig. ^a). This figure 
shows data from a Monte Carlo (MC) simulation of a two-dimensional ki- 
netic Ising ferromagnet below its critical temperature, which is driven by a 
sinusoidally oscillating field. However, the shape of the loops shown is quite 
general. For simplicity and concreteness, in this paper we use magnetic lan- 
guage, designating the oscillating force the "field" and the system response 
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the "magnetization" . These quantities can easily be re-interpreted when one 
discusses other types of systems, such as those mentioned above. 

Hysteresis was first systematically investigated in the late 19th century by 
engineers and physicists primarily concerned with the development of elec- 
tric motors and transformers The area of the hysteresis loop is 
proportional to the magnetic energy loss during one field period, as was first 
pointed out by Warburg [||. Its dependences on the frequency and ampli- 
tude of the applied field have therefore been studied intensively ever since, 
but open questions still remain. In particular, the question of whether the 
low-frequency behavior for ultrathin films of highly anisotropic materials is 
asymptotically a power law or logarithmic is still under investigation, both 
experimentally §,@,U|,|l§ and theoretically [|l]jl|,|l|,0jl5[|l|,|l3 . 



A different aspect of hysteresis in bistable systems, which is the main 
topic of the present paper, occurs at higher driving frequencies. When the 
frequency becomes sufficiently high, the symmetry of the hysteresis loop, 
which is apparent in Fig. 0(a), is broken. Instead of oscillating between its 
two stable values with a phase lag relative to the field, the magnetization 
oscillates around one or the other of its zero-field stable values. A series of 
such asymmetric loops is shown in Fig. [^(b). This symmetry breaking has 
become a topic of vigorous research during the last two decades. It was first 
reported by Tome and de Oliveira jist , who observed it during numerical 
solutions of a mean-field equation of motion for a ferromagnet in an oscillating 
field. 

Subsequently, symmetry breaking has been observed in numerous MC 
simulations of kinetic Isi ng sys tems, p^ , p^ , pO 21 2^ , ^ , p^ , as well as in fur- 



ther mean-field studies |20| , p2| , P3y25f| . It may also have been experimentally 
observed in ultrathin films of Co on Cu(OOl) Reviews of the field as 

it stood in 1994 and as it stands today can be found in and re- 
spectively. There now appears to be a consensus that the symmetry breaking 
corresponds to a genuine second-order, nonequilibrium phase transition. As- 
sociated with the transition is a divergent time scale (critical slowing-down 
[ p2| and, for spatially extended systems, a divergent correlation length ||l7| , p^ . 
Although estimates of the critical exponents have recently become available 
from finite-size scaling analyses of MC data for a two-dimensional Ising sys- 
tem in a sinusoidally oscillating field (lj,|2^ , their accuracy is not yet sufficient 
to decide whether the dynamic phase transition in this system belongs to a 
previously known universality class or represents a new one. Neither are we 
aware of theoretical arguments that can resolve the issue or help to identify 
the correct "critical clusters" . 

The main purpose of this paper is to present preliminary results for two- 
dimensional kinetic Ising systems in square- wave oscillating fields. These may 
help point the way towards answering some of the remaining questions con- 
cerning this intriguing nonequilibrium critical phenomenon. The use of a 
square-wave field both tests the universality of the dynamic phase transition 
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and significantly increases the computational speed. The structure of the re- 
mainder of the paper is as follows. The kinetic model system and the relevant 
quantities, including the dynamic order parameter, are defined in Sect. ^. A 
brief primer on metastable decay in spatially extended systems is given in 
Sect. ^ Numerical results are presented in Sect. ^, and a discussion and some 
suggestions for further research are given in Sect. |. 



2 Model and Relevant Quantities 

Most numerical simulations of the dynamic phase transition (hereafter abbre- 
viated DPT) in spatially extended systems have been performed on nearest- 
neighbor kinetic Ising ferromagnets on hypercubic lattices with periodic bound- 
ary conditions. These models are defined by the Hamiltonian 

n = -JJ2s.,S,-H{t)Y,S^, (1) 

where Si — ±1 is the state of the ith spin, ) runs over all nearest-neighbor 
pairs, J > is the ferromagnetic interaction, J^i runs over all L'^ lattice sites, 
and H(t) is an oscillating, spatially uniform applied field. The magnetization 
per site is 

m{t)^L-''Y,s,{t) . (2) 

i=l 

The temperature T is fixed below its zero- field critical value Tc, so that the 
magnetization for H = has two degenerate spontaneous equilibrium values, 
±TOsp(T). For nonzero fields the equilibrium magnetization has the same sign 
as H, while for not too strong H (see Sect. ^ for quantitative statements) 
the opposite magnetization direction becomes metastable and decays slowly 
towards equilibrium with time. 

The dynamic used here, as well as in and |24|, is the Glauber single- 



spin- flip MC algorithm with updates at randomly chosen sites 1 28 . The time 
unit is one MC step per spin (MCSS). Each attempted spin flip from Si to 
—Si is accepted with probability 

Here AEi is the energy change that would result from accepting the spin flip, 
and P = l/kBT where fce is Boltzmann's constant. For the largest system {L — 
512 square lattice) we employed a massively parallel implementation of this 
algorithm | |29| , p0||3lt|3^ . Other dynamics that have been used in MC studies of 
the DPT are Glauber or Metropolis |^ with updates at sequentially selected 
sites [ 20 , 2l],p2|,p3|p^,|27[| . Although the choice of update scheme can lead to 



subtle differences in the dynamics [B3| and we prefer random site selection 
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as the more physical scheme, we do not beheve it affects universal aspects of 
the DPT. 

The dynamic order parameter is the period-averaged magnetization JTsf , 
^ ' m(t) dt , (4) 



2t,/ 

where ti/2 is the half-period of the oscillating field, and the beginning of the 
period is chosen at a time when H{t) changes sign. Analogously we also define 
the local order parameter 

^ ' s,{t) dt , (5) 



which is the period- averaged spin at site i. For slowly varying fields the prob- 
ability distribution of Q is sharply peaked at zero [Ij . We shall refer to this 
as the dynamically disordered phase. For rapidly oscillating fields the Q distri- 
bution becomes bimodal with two sharp peaks near ±msp(T), corresponding 
to the broken symmetry of the hysteresis loops We shall refer to this 
as the dynamically ordered phase. Near the DPT we use finite-size scaling 
analysis of MC data to estimate the critical exponents that characterize the 
transition. 

Previous studies of the DPT have used an applied field which varies si- 
nusoidally in time. While sinusoidal or linear saw-tooth fields are the most 
common in experiments and are necessary to obtain a vanishing loop area 
in the low- frequency limit the wave form of the field should not af- 

fect universal aspects of the DPT. This should be so because the transition 
essentially depends on the competition between two time scales: the half- 
period ti/2 of the applied field, and the average time it takes the system to 
leave the metastable region near one of its two degenerate zero-field equi- 
libria when a field of magnitude Hq and sign opposite to the magnetization 
is applied. This metastable lifetime, {t{Hq,T)), is usually estimated as the 
average first-passage time to zero magnetization. In the present paper we 
use a square-wave field of amplitude Hq. This has significant computational 
advantages over the sinusoidal field variation since we can use two look-up 
tables to determine the acceptance probabilities: one for H = +Hq and one 
for H = -Hq. 

In terms of the dimensionless half-period, 

= ti/2/(r(ffo,T)) , (6) 

the DPT should occur at a critical value 0c of order unity. Although can 
be changed by varying either ti/2, Hq, or T, in a first approximation we ex- 
pect 0c to be independent of Hq and T. This expectation is confirmed by 
simulations carried out at several Hq and T for different system sizes. In par- 
ticular. Fig. |^(a) shows the average norm of the order parameter {\Q\) vs. 
for various field amplitudes and the corresponding metastable lifetimes. For 
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weaker fields (longer lifetimes) the transition is apparent at some 0^ ^ 1, 
while it clearly disappears (no dynamically ordered phase exists) for suffi- 
ciently strong fields (small lifetimes). Figure ||(b) shows Oc vs. {t{Hq,T)) 
where 0c was determined approximately as the location of the peak in the 
fluctuations of Q in a L = 64 system. Note that 0c for an infinite system is 
slightly different from this estimate due to finite-size effects. Also, these effects 
are expected to be smaller for strong fields (small lifetimes) where the spins 
become uncorrelated. The deviations from unity for very small {t{Hq,T)) 
and the complete disappearance of the DPT are discussed further in Sect. 

In many studies of the DPT the transition has been approached by chang- 
ing Ho or T fl^,^ 21,2^. While the above discussion indicates that this is 
correct in principle, we show in Sect. |^ that {t{Hq, T)) depends strongly and 
nonlinearly on its arguments. We therefore prefer changing at constant 
Hq and T as this in practice gives more precise control over the 

distance from the transition. 



3 Decay of Metastable Phases 

In mean-field studies of the DPT Hq must be larger than the temperature- 
dependent spinodal field beyond which the metastable free-energy minimum 
disappears For weaker fields the hysteresis loops remain asymmetric, 
even for the lowest frequencies. It is much less appreciated that bounds on 
the fields and temperatures for which a DPT can occur also exist for spa- 
tially extended systems. These bounds are readily obtained from classical 
nucleation theory and the Kolmogorov-Johnson-Mehl-Avrami (KJMA) 
theory of metastable decay 0,||,|3[||,||,|3 by comparing four character- 
istic lengths. These are the lattice constant (here defined as unity), the size 
i?c(|^^|, T) of a randomly nucleated critical droplet of stable phase, the typical 
distance Rq{\H\,T) between individual droplets, and the system size L. 

The critical radius is determined by the competition between the sur- 
face free energy of the droplet, oc (t(T)R'^~^ where (j{T) is the surface ten- 
sion between the two phases, and the bulk free energy, oc —\H\R'''. As a 
result, Rc ex. a{T)/\H\. The nucleation rate per unit volume is I{H,T) oc 
exp [— S'(r)/|i7|''~^] , where S{T) is the field-independent part of the free 
energy of a critical droplet divided by /cbT, and preexponential powers of 
have been suppressed. 

The classical KJMA theory describes the metastable decay as homoge- 
neous nucleation of droplets at random times and positions with nucleation 
rate /, followed by deterministic growth of these droplets with constant radial 
velocity v oc Assuming that the droplets overlap freely when they meet, 
the time dependent magnetization is given by the well-known "Avrami's 
Law", 



m{t) « rrigp < 2 exp 
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2 exp 



d+i 



- 1 



(7) 



where Qd is a proportionality constant such that the volume of a droplet of 
radius R equals fidR"^- The argument of the exponential is the "extended vol- 
ume" [ p6| , i.e., the total volume fraction of stable- phase droplets, uncorrected 
for overlaps. Solving (|^) for the time at which m — gives the lifetime. 



{r{\HlT)) 



l)ln2 



(8) 



which depends on \H\ and T through v and /, but is independent of L. The 
characteristic length i?o is obtained from v{\H\,T) and {t(\H\,T)) as 



Ro{\H\,T) ^ v{\HlT) (r(|i/|,r)} cxexp 



(9) 



where a preexponential power oi\H\ again has been suppressed ]3^ , p5{ . 

The regime in which a DPT can occur is that in which a large number of 
droplets contribute to the decay of the metastable phase: 



K i?c < ^0 < 



(10) 



This is known as the multidroplet (MD) regime |^,^. It is limited on the 
weak- field/small-system side by the dynamic spinodal (DSP) field , 



1 ^{T) 
d+l InL 



1/(^-1) 



(11) 



which corresponds to Rq w L. For \H\ < i?DSP almost always only a single 
droplet contributes to the magnetization switching. A DPT does not occur 
in this regime, but stochastic resonance is observed at low frequencies pSf . 
On the strong-field side the MD regime is limited by an L-independent cross- 
over field approximately given by 2Rc « 1. By a somewhat confusing term this 
crossover is often referred to as the "mean-field spinodal" i?MFSP (T) p^ , |35{ . 
In the strong-field (SF) regime beyond this limit the spins become increas- 
ingly uncorrelated as \H\ increases. This is the region of small {t{Hq,T)) 
(large Hq) where 0c rapidly approaches zero and the transition disappears 
(Fig. |2|). In computational studies of the DPT it is essential to ensure that 
HnSp{T,L) < Ha < HupspiT) for all values of L and T used. 



4 Results 

We performed extensive simulations on square lattices with L between 64 
and 512 at T — 0.8Tc and Ho = 0.3J. Typical run lengths near the transition 
range from 2x10** to 10^ full periods, corresponding to 2.8x10^ — 1.4x10*^ 
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MCSS. The system was initialized with all spins up and and the square- 
wave external field started with the half-period in which H = —Hq. After 
some relaxation the system magnetization reaches a limit cycle (except for 
thermal fluctuations), i.e., Q is stationary. We discarded the first 500 periods 
of the time series to exclude transients from the stationary-state averages. 

For small half-periods (6* ^ 0c) the magnetization does not have time 
to switch, resulting in ~ mgp, while for large half-periods (0 » 0c) it 
switches every half-period and Q « as can be seen from the time series in 
Fig. Ij. The transition between the high- and low-frequency regimes is singu- 
lar, characterized by large fluctuations in Q. To illustrate the spatial aspects 
of the transition we also show configurations of the local order parameter 
Qi in Fig. ^. Below 0c [Fig. ^a)] the majority of spins spend most of their 
time in the -1-1 state, i.e., in the metastable phase during the first half-period, 
and in the stable equilibrium phase during the second half-period (except for 
equilibrium fluctuations). Thus, most of the Qi w +1. Droplets of Si — — 1 
that nucleate during the negative half-period and then decay back to -1-1 
during the positive half-period show up as roughly circular gray spots in the 
figure. Since the spins near the center of such a droplet become negative 
first and revert to positive last, these spots appear darkest in the middle. 
Above 0c [Figs. §(c,d)] the system follows the field in every half-period (with 
some phase lag) and Qi~0 at all sites i. Near 0c [Fig. |4|(b)] there are large 
clusters of both Qi « -f 1 and Qifn—l separated by wide "interfaces" where 
Qi wO. Also, for not too large lattices one often observes the full reversal of 
an ordered configuration {Qi} —>■ —{Qi}, typical of finite, spatially extended 
systems undergoing symmetry breaking. 

For finite systems in the dynamically ordered phase the distribution of Q 
becomes bimodal. Thus, to capture symmetry breaking, one has to measure 
the average norm of Q as the order parameter, i.e., (|Q|). Figure ^(a) clearly 
shows that this order parameter is of order unity for < 0c and vanishes 
for > 0c, except for finite-size effects. 

To characterize and quantify this transition in terms of critical exponents 
we employ the well-known technique of finite-size scaling ||28| , ^ . The quantity 
analogous to the susceptibility is the scaled variance of the dynamic order 
parameter. 



^L-i'((Q')L-(|Q|)i) . (12) 

For finite systems has a characteristic peak near 0c [see Fig. 1(b)] which 
increases in height with increasing L, while no finite-size effects can be ob- 
served for <^ 0c and :s> 0c- This implies the existence of a divergent 
length scale, possibly the correlation length which governs the long-distance 
behavior of the local order-parameter correlations (QiQj). Note that the lo- 
cation of the maximum in Xl shifts with L. This also contains important 
information about the critical exponents. 
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To estimate the value of 0c at which the transition occurs in an infinite 
system we use the intersection of the fourth-order cumulant ratios p^ , |39t , 

for several system sizes as shown in Fig. ||(c). For the largest system {L = 
512) the error bars on Ul were too large to use it to obtain estimates for 
the crossing. Our estimate for the dimensionless critical half-period is = 
0.913 ± 0.003 with a fixed-point value U* = 0.615 ± 0.005 for the cumulant 
ratio. 

For our model the quantity analogous to the reduced temperature in equi- 
librium systems (i.e., the distance from the critical point) is 

Finite-size scaling theory provides simple relations for the order 
and its scaled variance Xl for finite systems in the critical regime 

{\Q\)l = L-^'/''T±{eL^''') , (15a) 
Xl ^ L^/^g^ieL^/") , (15b) 

where J-'± and g± are scaling functions and the -f (— ) index refers to 6* > ©c 
{0 < 0c)- Then at 0c {9 — 0) we have the finite-size behavior of the above 
two observables: 

{\Q\)l oc L^^/^ , (16a) 
Xl oc L'^/'^ . (16b) 

Using {\Q\)l and Xl for several system sizes a t 0c (estimated as the value 



of where the cumulants cross) we employ (16a) and (16b) to find the 
exponent ratios fS/v and "//ly through weighted linear least-squares fitting to 
the logarithmic data. Further, the shift in the location of the peak in X^ for 
these finite systems, 0c{L), is used to estimate v in the same way ||2^ , ^ : 

\0c{L)-0c\^L-^''' . (17) 

Our estimates for these exponents are (3/v=Q.122 ± 0.005, ^jv=\.ll ± 0.05, 
and z/=1.0 ± 0.15. The largest uncertainty occurs in v due to the large rel- 
ative errors in |0c(i) — 0c\- We note that these numbers are very close 
(within the one-standard-deviation error bars) to the critical exponents of 
the two-dimensional Ising universality class /3/^=l/8=0.125, 7/i^=7/4=1.75, 
and v=\. However, as discussed in Sect. ||, we do not consider this conclusive 
proof of the universality class for this noncquilibrium phase transition. 



To obtain a more general picture of how well the scaling relations in (15a) 
and ( p^ are obeyed, we plot (\Q\)lL^I'' [Fig. |(a)] and XlL-^^" [Fig.|§(b)] 



vs. eP^ ra. In these fi gures we used the exponents of the two-dimensional 
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Ising universality class since they are within one standard deviation of our 
estimates and analogous data-collapse plots using our numerical exponent es- 
timates do not look perceptibly different. The plots show excellent agreement 
with the scaling assumption and graphically define the scaling functions J-± 
and Q± . The asymptotic behaviors of these functions have to be simple power 
laws so that the true critical behavior is restored when the limit L — > oo is 
taken in (f5a) and (15b). Some deviation can be observed in Fig. ||(b): below 
0c the data points for the smaller systems systematically start to peel off 
earlier from the straight line representing the asymptotic behavior of for 
large argument. 

In the absence of theoretical arguments why this DPT should belong to 
the Ising universality class we also tested the data collapse with the expo- 
nents for random percolation, which are relatively close to the corresponding 



Ising ones gl]: /3/i/=5/48«0.104, 7/i/=43/24wl.79, and i/=4/3«1.33. It ap- 
pears reasonable to consider this universality class in particular, since the 
metastable decay process described in Sect, y produces transient spanning 
clusters that belong to the random-percolation class as m{t) passes through 
H- and T-dependent percolation thresholds for the two phases near m — Q 

. Data-collapse plots using the random-percolation exponents are shown 
in Figs. ^(c,d). They are clearly inferior to the ones with the Ising exponents 
[Figs. ^(a,b)]. In particular, for > 0c (the dynamically disordered phase) 
the data collapse gets progressively worse for larger systems. Comparing the 
scaling plots in Figs. |(a,b) with Figs. ^(c,d) we conclude that the exponent 
v for the DPT is significantly different from that of the random percolation 
universality class and that it is closer to the Ising value. This conclusion is 
supported by direct comparison of the random-percolation and Ising expo- 
nents with our numerical estimates for the DPT. Our estimates for p/iy and 
u lie more than two standard deviations away from the percolation values, 
but less than one standard deviation away from the Ising values, while our 
estimate for 7/1/ lies midway between and less than one standard deviation 
away from the exact values for both classes. 

The exponents previously obtained with a sinusoidally oscillating field 
[ p"7||3^ lie within two standard deviations of our values. As it was difficult to 
determine the uncertainties in the estimates obtained in that study, we believe 
the results are consistent and indicate that the DPT in systems driven by 
sinusoidal and square-wave fields belong to the same universality class, as 
expected. 



5 Summary and Discussion 

In this paper we have presented preliminary results from a large-scale Monte 
Carlo study of the dynamic phase transition (DPT) in a two-dimensional 
kinetic Ising ferromagnet driven by a square-wave oscillating field. Our re- 
sults are consistent with those of previous studies of the same model in a 



10 



sinusoidally varying field [p?7|j24[ | . They indicate that for field amplitudes Hq 
such that the metastable magnetized phase decays to equilibrium via the 
multidroplet (MD) mechanism described in Sect. the system undergoes a 
continuous DPT when the half-period of the field, ti/2, approximately equals 
the metastable lifetime, {t{Hq, T)). Thus the critical value 0c of the dimen- 
sionless half-period defined in (^) is near unity. As O is increased past 0c 
the order parameter which is the expectation value of the norm of the 

period-averaged magnetization, vanishes in a singular fashion, as shown in 
Fig. 1(a). 

The strong and systematic finite-size effects in the order parameter and its 
scaled variance indicate that there is a divergent correlation length associated 
with the transition. Using standard finite-size scaling techniques borrowed 
from the theory of equilibrium phase transitions we therefore estimated 0c 
and the critical exponents (3, 7, and 1^ from data for system sizes between L = 
64 and 512. The resulting estimates are /3/j/=0.122 ± 0.005, -f/v=1.77 ± 0.05, 
and i/=1.0 ± 0.15, which agree to within the statistical errors with those 
previously obtained with a sinusoidally oscillating field ||l^,^ . This is strong 
evidence that the shape of the wave does not effect the universal aspects of 
the DPT. 

We also note that our estimates agree well within one standard devia- 
tion with the exact values for the two-dimensional equilibrium Ising model: 
/3/jy=l/8=0.125, 7/1^=7/4=1.75, and v=l, and that they satisfy the expo- 
nent relation 



2(/3/t/) + (7/1^) = 2.01 ±0.05 «d. (18) 

The fixed-point value of the fourth-order cumulant ratio, U* =0.615 ± 0.005, 
is also close to that of the Ising model, C/* = 0.610 690 1(5) [Esl. As discussed 



in 1 17 we find the estimated values for the exponents and U* conclusive 
evidence that the DPT corresponds to a nontrivial fixed point, so that ( p^ 
indeed represents a hyperscaling relation. 

Since we know of no convincing theoretical argument that this nonequi- 
librium phase transition should be in any particular equilibrium universality 
class, we also compare the exponent estimates with those of the universality 
class of random percolation. For this class the values ^l[ /3/z^=5/48«0.104 
and i^=4/3«1.33 are significantly different from our exponent estimates, while 
7/j/=43/24wl.79 lies within the uncertainty in our estimate. Taken together, 
the three estimates and the low quality of the data collapse plots shown in 
Figs. H(c,d) constitute strong evidence that the DPT is not in the universality 
class of random percolation. However, we emphasize that we do not yet con- 
sider the question of the universality class as settled. That will require still 
more accurate numerical exponent estimates as well as theoretical arguments. 

Outside the MD regime the phase transition disappears. For weaker fields 
or smaller systems the metastable phase decays stochastically through a single 
droplet, and no long-range correlations evolve 138]. In the other extreme, as 
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Hq is increased into the strong-field regime, 0c goes to zero and vanishes at 
a sharply defined field (lifetime), as shown in Fig. ^(b). 

In conclusion, the dynamic phase transition observed in spatially extended 
kinetic Ising systems driven by oscillating fields is a fascinating nonequilib- 
rium critical phenomenon with all the hallmarks of a thermal phase transition 
corresponding to a nontrivial fixed point. Although the present study con- 
firms that the particular wave form of the oscillating field does not change 
the universality class, and our numerical estimates for the critical exponents 
and cumulant ratio are consistent with the universality class of the two- 
dimensional Ising model in equilibrium, many questions remain. These in- 
clude the question of universality under changes in T and Hq , and the details 
of the disappearance of the transition as Hq is increased. No procedure to 
identify critical correlated clusters from the connected clusters seen in con- 
figuration "snapshots" such as Fig. ^(b) is yet known, and theoretical ideas 
are largely missing. Thus, there is still much to be done! 
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Fig. 1. Series of hysteresis loops obtained from a Monto Carlo (MC) simulation of 
a two-dimensional kinetic Ising ferromagnet at r = 0.8 Tc in a sinusoidally varying 
field. The dimensionless magnetization m is shown vs. the field H in units of the 
nearest-neighbor interaction energy J. (a) : A series of regular symmetric hysteresis 
loops, typical of a system driven at a moderately low frequency, (b): A series of 
asymmetric hysteresis loops, typical of a system driven at a high frequency 
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Fig. 2. (a) Order parameter vs. O for L = 64 square system at r=0.8Tc for 
several values of the field amplitude Hq. The corresponding metastable lifetimes, 
{t{Ho,T)) given in MCSS, are shown in the legends, (b) The approximate dimen- 
sionless critical half-period 0c vs. the metastable lifetime {t{Ho,T)) in MCSS 
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Fig. 3. Time scries of the order pa- 
rameter Q at T = 0.8Tc and Ho = 0.3 J 
for L = 128. Horizontal trace near 
Q = +1: 6> = 0.27 < 0c (dynami- 
cally ordered phase). Strongly fluctu- 
ating trace: 6 = 0.98 « Oc (near the 
DPT). Horizontal trace near Q = 0: 
= 2.1 > ©c (dynamically disordered 
phase) 




Fig. 4. Configurations of the local order parameter Qi at T = 0.8Tc and Ho = 0.3 J 
for L= 128. (a) = 0.27 < 6»c (dynamically ordered phase) . (b) 6» = 0.98w6>c (near 
the DPT), (c) 6> = 2.7 > 6>c (dynamically disordered phase), (d) 6> = 13.4 > 6>c 
(deeper in the disordered phase). On the gray-scale —1 (+1) corresponds to black 
(white) 
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Fig. 5. Finite-size behavior of our ob- 
servables at T = 0.8Tc and Ho = 0.3 J for 
various system sizes, (a) The order pa- 
rameter (IQI)l. (b) The scaled variance 
Xl as defined in (y). (c) The fourth- 
order cumulant ratio as defined in 
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Fig. 6. Finite-size scaling plots at T = O.STc and Ho = 0.3J using the two- 
dimensional Ising exponents: (a) For the order parameter (|Q|)l. (b) For the scaled 
variance Xl- Finite-size scaling plots at the same temperature and field using the 
random-percolation exponents: (c) For the order parameter (|Q|)l. (d) For the 
scaled variance Xl - The straight lines in all graphs represe nt th e asy mpto tic large- 
argument behaviors of the functions T± and C}± given in (15a) and (|l5b|), respec- 



tively. The value for the (infinite-system) 0c, used in all four graphs, was obtained 
from the cumulant crossing [(|l3|)] Fig. ^(c) 



